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I. INTRODUCTION 



Hyperon radiative decays constitute an interesting class of processes that has been stud- 
ied intensively, both theoretically and experimentally, for the past 20 years. However, despite 
these efforts we still lack a global understanding of these processes [0. 

Some of these decays may be sensitive to the s ^ d'j transition, which tests the Stan- 
dard Model at the one-loop level and is similar to the b ^ s'y transition recently seen at 
Cornell which has generated a large amount of interest . Therefore, it would be extremely 
interesting to obtain any experimental information on it. In this paper we re-examine the 
assertion ^ that the decay can provide a window to the s d'y transition. In 

the following we briefly review the arguments that led to this conclusion. 

There are two distinct contributions to hyperon radiative decays, namely short- and 
long-distance contributions. In short-distance processes, the photon is emitted at a close 
distance (typically x ~ 1/Mw) from the weak interaction vertex whereas in long-distance 
processes [x ~ confinement radius) one can separate the weak from the electromagnetic 
vertices. 

At the quark-level, there are in general three types of contribution to hyperon radiative 
decays, as shown in Fig. 1. The first diagram (Fig. la) corresponds to a W-exchange and 
is a long-distance process, whereas Fig. lb (single-quark transition s — >• dj) and Fig. Ic 
(penguin diagram) are short-distance processes. 

It is known that the s ^ d'j transition cannot be the dominant one in all hyperon decays 
0. This is to be expected since W-exchange is the main term in the AS" = 1 effective 
Hamiltonian. However, there are two hyperon radiative decays, namely fl~ — >■ S~7 and 
S~ —>■ S~7, where the quark content of the initial and final state baryons does not allow a 
W-exchange contribution. Therefore, only the s —>■ d'j and penguin diagrams contribute at 
the quark level to the decay. 

Let us concentrate first on ^"7. It has been shown that the penguin contribu- 

tion for this process is negligible. In Ref. [Q, the s — > 0^7 contribution to the branching ratio 
was calculated to be BW^'^^{'Er~ — > S~7) = 1.8 x 10~^. However, although the W-exchange 
process does not contribute directly to 'Er S^7, it does contribute to — > tt^A, where 
the final state can then re-scatter as tt^A ^"7. This long-distance contribution has 
been estimated as BR^^{E~ ^"7) = 1.7 x 10~^ 0], close to the experimental result 
BRE^p(E- S-7) = (1.27 ± 0.23) X 10"^ §. 

This leaves us with the process —>■ S~7 as the only candidate for probing the s d'j 
transition in hyperon radiative decays. Let us summarize what is known about this process 
so far. There is an experimental upper limit determined recently 0: 

BR^^'P-iQ- H-7) < 4.6 X 10"^ (1.1) 

This is well above the unitarity limit : 

The penguin contribution has been estimated as |0 : 
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Kogan and Shifman have also calculated the long-distance contribution from the dominant 
E^n' intermediate state as : 

BR^^{Q- E--f) = (1-1.5) X 10~^ (1.4) 

The contribution from the single-quark s ^ d'y transition has been re-evaluated recently 
ig yielding : 

BR'^'^^Q- E--f) = 8.0 X 10"^ (1.5) 

This result corrects a previous calculation [Q, where it was claimed that the single quark 
transition should be the dominant contribution to the Q~ radiative decay. Both calculations 
10|,|3[| were based on estimating the hadronic matrix element of the effective operator de- 



scribing the single-quark transition, using SU(6) quark- model wave functions and assuming 
unit overlap. 

Due to the importance of this process, we believe that a more reliable estimate of the 
single-quark s ^ d'j contribution to the radiative decay H~7 is in order. This is the 

purpose of the present work, where we perform a calculation of the same matrix element 
using the well known technique of QCD sum rules |]TT] . 

The QCD sum rule approach was successfully applied to the E"*" — > p7 decay by Balitisky 
and collaborators |T^, who calculated in this way the dominant W-exchange. QCD sum- 



rules were also employed before for the calculation of the single-quark contribution [|T^ to 
the S+ and S~ radiative decays. Although the single quark transition is not the dominant 
contribution in these decays, the calculation exposed the role of non-perturbative corrections 
which alter the naive picture of the single-quark transition in a major way: it showed that it 
is not a good approximation to assume a unity overlap for the quark-model wave functions. 

This paper is organized as follows. Section II discusses the effective operator for the 
s — > (i7 transition. In Section III we derive the sum rules and in Section IV we present the 
results and conclusions. 



II. THE EFFECTIVE OPERATOR FOR S ^ D-f TRANSITION 

We are interested in computing the contribution of the quark level s ^ d'y transition to 
the radiative decay Q H7. This transition does not occur at tree-level in the Standard 
Model. Hence, it could in principle provide a good test to the Standard Model at the 1-loop 
level, being sensitive also to new Physics. 

The effective hamiltonian (with the heavy quarks c, b and t as well as electroweak gauge 
bosons integrated out ) that describes \AS = 1| transitions is given by : 

^eff"^" = ^ cMOM , (2.1) 

where Gp is the Fermi constant, and we use the notation Aj to denote the following product 
of elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, Aj = VsiV*^ . The Wilson 
coefficients Ck can be computed perturbatively and {Ok} is a complete set of operators 
written as: 
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q=u,d,s 
q=u,d,s 

O, = (d^^^PLS^) E iQPl'^PRQfs) 

q=u,d,s 



q=u,d,s 



le 



167r2 
IGvr 



ms{daaf,^PRSa)F^'' 
^ms{d^a^un^PRSf,)G''^\ 



(2.2) 



where P l = ^(1 =F 75) , o'au = ^[Tu, lulf^] F^'^ and G"-^'^ are the electromagnetic and color 

R 

field tensors respectively. 

We are interested in the Wilson coefficient Ci{fi = rris), which controls the short distance 
contribution to the process s —>■ d'j. 

The Wilson coefficients are first computed at the scale /i = M\y in zeroth order in QCD. 
QCD corrections are included by evolving these coefficients down to /i = using the renor- 
malization group equations. This evolution is accomplished in three steps. Schematically 
one has : 



Cfc(/i = iris) = Uli{ms,mc) U^^{mc,mb) U^^{mb, Mw) Cm{Mw)- 



(2.3) 



The evolution matrix with / active quark ffavors in the leading logarithmic approximation 
can be written as : 



f/^(mi,m2) = VTV-^ 



(2.4) 



where is a matrix that diagonalizes the transpose of the 8x8 anomalous dimension matrix 
7 : 

Id = V~'^V , 
and T is a diagonal matrix with elements given by : 

as{m2) 



(2.5) 



a^(mi) 



(2.6) 



where (3q = 11 — 2//3 and in this approximation the running of the strong coupling constant 
is described by : 



4 



§0 
2 



^a^(mi)ln(mi/m2) 



(2.7) 



We now turn to a discussion of the initial conditions. Due to our choice of factoring out 
Xu in eq.( p.l|) , we have: 



C2(Mi 



w 



Cl,3-6(M, 



. 



(2.8) 



The values for cj and cs at /i = Mw are given in terms of functions arising from one-loop 
diagrams with an external photon and a gluon respectively |jl5|,|l6|. These function F2{xi) 
and D{xi), where Xi = mf/M^r and is the mass of the quark running in the loop, are 
given by: 



F2{x) = Q 



x^ — — 2x ^ 3x^ \n{x) 



D(x) 



4(x - 1, 
3x^ hi{x) 
2{x-iy ' 
x^ — — 2x 



\3 



4(x 



+ 



2{x 



ln(x) 
2{x-iy 



+ 



2x^ 



4(x 



(2.9) 
(2.10) 



where Q is the charge of the internal quark. 

In table 1 we show the contributions to F2 from the different internal quarks and its 
product with the relevant CKM matrix elements. We use the following values for the quark 



masses and A,- 



nit 



175 GeV, rric = 1.5 GeV, m„ = 5 MeV, A„ = 0.21, A, 



-0.21 



and At = —(1.2 — 7.2) x 10""^, where the CKM elements were obtained from the standard 
parametrization of the Particle Data Book neglecting CP violation, that is taking ^13 = 0. 

From this table we see that, as opposed to the b ^ sj case, the top quark contribution 
is not the dominant one for s —>■ d'-f and we must take into account both t— and c— quarks 
in computing the initial values for C7 and cg, which are given by : 




hF,(xt) + ^F2{x^ 
^D(Xt) + ^D(x, 




(2.11) 
(2.12) 



In order to compute the evolution matrix we use the anomalous dimension matrix given in 
ref. fl^ and as^M^] 



0.124. However, from eq. ( |2.7|) one would get a value for as{ms) > 1, 
since for consistency with the leading-log anomalous dimension matrix we have to compute 
the running of in one-loop only. In this case we use as{ms) = 1. We find: 



crirris) = -0.50c2{Mw) + Q.27cj{Mw) + Q.l3cg{Mw). 
The numerical values for the initial conditions for c-j and Cg are: 

Cj{Mw) = (2.1 - 7.7) X 10"^ , CsiMw) = (1-4 - 4.2) x 10^^ 



(2.13) 



(2.14) 
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and from Eg. ( |2.13| ) we see that ct^ttIs) is insensitive to the the values of ct^Mw) and cs{Mw), 
being dominated by C2{Mw). This result is in contrast to what happens in the b ^ sj 
case. It can be understood since in the s ^ d'y case we have C2{Mw) oc A„ and cj^s{M\y) is 
effectively proportional to with At <^ in the b ^ case, one has C2{M]y) oc iHcKtl — 

\VbtVZ\ (X C7,8iMw) . 

Unfortunately, Physics beyond the Standard Model would contribute mostly to cj^siMw) 
and from our calculation we see that unless the modification implies in an increase of 
Ct^s{Mw) by roughly four orders of magnitude, no significant contribution will be made 
to the s — > (i7 process. 

Using equations (|2.2|), ( p.8D and ( |2.13| ) we arrive at the main result of this section: 



^eff.'" = Crim,)Or = (0.28) ^m,(4cr^,PHS„)F'^^ . (2.15) 



III. QCD SUM RULES 

The subject of this section is the evaluation of the s — c/7 contribution to the four 
multipole amplitudes that appear in the transition matrix related to the radiative decay 
S~7, using a QCD sum rule approach. 

A. Formalism 

We start with the three-point correlation function 

l^,{Q^,Qf) = rfWi/e^«/^e-«'-(0|T[r/=(i/)K;/"(0)rJ^(x)]|0) , (3.1) 

where 'H'^ejf^^ given by Eg. ( |2.15| ). 

As usual, our goal is to make a match between the two representations of the correlation 
function ( p.lD at a certain region ~ 1 GeV^: the operator product expansion (OPE) in 
powers of and the phenomenological representation of the dispersion integrals. The basic 
idea, supported by ample successful applications, is that taking into account only the first 
few terms in the Qj expansion, complemented by rather simple assumptions for the higher 
mass contributions to the dispersion relation, will already provide a good estimate to the 
amplitudes of interest. 

The VL~ and S~ interpolating fields are given by []TB 



T]^ {x) = -eabc (^sl{x)C-ff,Sb{x)^ lbl^dc{x) , (3.2) 

(x) = ^abc [s^{x)Cj^Sb{x)) Sc{x) . (3.3) 

With these definitions and working at leading order in perturbation theory we obtain for 
the T product appearing in Eg. (|0|) 
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{trlS^Ay - xKC{S:Ay - x)fC^u] + 2^,C{S:Ay - x)fC^,Sl,{y - x)} , (3.4) 

where the coordinate-space quark propagator in the presence of the quark condensate takes 
the following form [T^JT^ 



Sl{x) = (0|T[g„(x)g,(0)]|0) = - - ^(gg) + '^rn,{qq)^ + ■ ■ ■ . (3.5) 



The amplitude ^^{Qi,Qf) includes a lot of different non-trivial Lorentz structures. For 



each one of these structures, Ilk{Qf, Q'i), we can write a double dispersion representation of 



the form 



n.(QlQ}) ^I,f,s,f,s,^-^I0^^ (3.6) 

where the ellipsis represents subtractions polynomials in Q? and , which will vanish under 



the double Borel transform pO[, which is a straightforward generalization of that used in 



Ref. [|TT|. Applying the double Borel transform to Eq.( p.6|) gives 

MMlMl) = - dsj d.2P.(5i,52)e-^^/^^'e-^^/^^ (3.7) 

TT^ JO Jo 

In the phenomenological side the various Lorentz structures can be obtained from the 
consideration of the Q and S contribution to the dispersion sum rule 

(o|r]1s(g;))(s(g/)7lK;/1^^(Q.))(^W.)l^^|o), (3.8) 

where the most general, gauge invariant form for the amplitude of the fl^ — > 2^7 decay is 
given by 

M{Q- ^ S-7) = {E{Qf)^\n:y;^\Q{Q,)) = ieGpu^^\QM(^i + ^^<^2m'"' -Yqn 

+ {hi + l,h){Q^.qg^' - QUn/M^]uf\Q::)e, , (3.9) 

where q = Q f — Qi, is the polarization of the photon, and ai, 02, 61 and 62 are the four 
amplitudes we want to evaluate. The other matrix elements contained in Eq.( |3.8|) are of the 
form 

(0|r|H(g^)) = AHn(=)(Q;), (3.10) 

MQ.)l^|0) = An<)(Q.), (3.11) 
where u^{p) is a Rarita-Schwinger spin- vector satisfying 

uf\p)u^:'\p) = - - ^7.7. - ^P,Pu + ^^^;~^^^" ) + Mn) , (3.12) 



u{p) is a Dirac spinor and Ah, Aq are the couplings of the currents with the respective 
hadronic states. 

Saturating the correlation function Eq. ( p.l| ) with Q and S intermediate states, and using 
Eqs.dp), (pD, (|3l0|) , (|3TTD and ( ^121) we get 



nf^")(g.,Q/) = 



Q2 + Ml 



,ai + 7502) 



+ (&1 +75&2)X 



3M, 



f7 



3M< 



+ Ml ■ 



(3.13) 



In this work we will derive sum rules for the Lorentz structures: 



(3.14) 
(3.15) 



a : ^(fli + 7502)0-0 
b : ^(61 + 75fe2)(Qi-ge^ - Qj.eg^) . 

This choice is based on the fact that the structure a gets perturbative contribution (which is 
very important in the case of double dispersion relation analysed here ||2^) as well as power 
corrections contributions, and the structure b is the only one that contributes solely to the 
amplitudes bi and 62- Furthermore, as we will show, both sum rules are very stable as a 
function of the Borel masses. 

For the continuum contribution we adopt the standard form of Ref . pO| , which completes 
our parametrization of the spectral density pk'- 



—pa[Si,S2, 



—2 

— eG'FAsAn5(si 



1 
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Ml)5{s2-Ml) 

+ ^(Sl-SH)^(s2-Sf,)pW(Si,S2) , 

= -ieG F>^B>^n5{si - Ml)5{s2 
+ 0{si - s-e)9{s2 - sn)pf\si, S2) 



Ml) 



(3.16) 



(3.17) 



where sh, are, respectively, the continuum threshold of the H and VL determined in the 
mass sum rules, and p^^\si, S2) is the free-quark spectral function. 

Let us now calculate the left-hand side of the sum rule Eq. (|3.7| ) using the OPE. We shall 
consider only the diagrams shown in Fig. 2. The perturbative contribution to the sum rule 
is shown in Fig. 2a and is obtained by using only the first term in the right-hand side of 
Eq.(|3.5D in Eq. (|3.4j). The resulting expressions for the two chosen structures are 



^(f)(Ml^M 



-3eGi70.28m, 



n(f)(M2,M2) 



2%6 (Mf + M|)3 

2%6 Im[TmI7 



f(l-75)^^a/37M^.e"g^ 



(3.18) 
(3.19) 



Noticing that Eq. (|3.7|) is the double Laplace transformation in 1/Mf, the free-quark 
spectral function, pf\si,S2), can be obtained from Eqs.( p.l8D and ( p.l9| ) by applying the 
inverse transformation. It gives 



8 



ipfe(Sl,S2) 



-3eG'ir0.28m, 



2107^6 

—ieGp0.28ms 3^, 
2%6 



SiS25(si - S2) 
sl6'{s2 - Si) . 



(3.20) 
(3.21) 



The next term in the OPE is determined by Figs. 2b, c and d. The contribution of Fig. 2d 
exactly cancells the contributions of Figs. 2b and c for the structure b. The result of the 
three diagrams for the structure a is 



-ms{ss} 



2H^ 



{Ml + M| 



2\2 



Ml + Ml In 



'Ml + M| 



;l-75)^^a/37M^.e"g^ (3.22) 



The last diagram which we take into account is shown in Fig.2e. It does not give any 
contribution to the structure b and the contribution to the structure a is 



-6^^0.28771, 

12^^2 



{ssfi{i - 75)^^^/37^^ 



(3.23) 



The gluon condensate is not taken into account because its contribution is suppressed relative 
to the operator mq{qq) (which has the same dimension) by an extra loop factor l/lGvr^. 

In general, when a real photon is emitted, one has also to consider diagrams such as that 
in Figure 3, involving long distances in the photon channel. However, this diagram gives a 
negligible contribution, as it is proportional to the part of the photon wave function at long 



distances, involving quarks of different flavors |T2|,|2T[. Therefore, we do not consider it here. 

Collecting all the obtained contributions, and transferring the continuum contribution 
to the OPE side, we arrive at the following representation to the amplitudes: 



3 e^^i/^^?e*^^/^2 



ai 



-a2 



+ 



27r2 

57r2 
16 



-0.28m, 



26 Jo Jo 



nisfa 



Ml Ml 



{Ml + M|)2 



Ml + Ml In 



ds2SiS2S{si - S2)e 
Ml + Ml^^ 1 



12'^ 



(3.24) 



1 Ml/Ml M^/M^ .g- rsn rl 

^1 = = 1^2 —^, 0.28m, rdsj ds2sl-^{6{s2 - s,))e-^^l^'h-^-l''^- 

l^T^ AhAq "'0 JO as 



(3.25) 



where \u = An Xh, a = — 47r (gg), and / = (ss) / (qq) 



B. Sum Rule Analysis 

The sum rules are sampled in the region of Borel mass, M^, which have been identified 



as the fiducial region for the baryon mass sum rules 

1.4 < Ml < 2.0 GeV^ 



for S , 



2.4 < M^ < 3.2 GeV^ 



forfi 



(3.26) 
(3.27) 
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The values used for and sh are extracted from the respective mass sum rules analysed in 
the same region given above. For S we get from Ref. P2[ — 1.6GeV^ and s= = 3.6GeV^. 



For Q, using the sum rules given in [^], we get Aq ~ 3.8GeV and ^ 7.5GeV . The value 
of the other parameters are H: a = 0.55 GeV^ m, = 150MeV, / = 0.8, Ms = 1.32GeV 



and Mn = 1.67GeV. 

In Fig. 4 we show the result obtained for function of for different values of 

M|. One can see that Oi varies very slowly with but not so slowly with M|. In the 
Borel region considered the amplitude Oi change less than 20%. The same behaviour can be 
observed in Fig. 5 where Oi is plotted as a function of M| for different values of M^. 

For the amplitude bi we get an even more stable result. As can be seen by Figs. 6 and 7 
bi varies slowly with M| and with Mf. In the Borel region considered bi changes less than 
10%. Using Mf = M| and M| = we get 

= -a2 ~ l.lOMeV, (3.28) 

6i = -62 ~ -0.58MeV. (3.29) 



IV. RESULTS AND CONCLUSIONS 



Given the values of the form factors 01,02,61,62, we can compute the decay width 
Tin- ^ S-7) i: 



where 



{Ml - Ml 



2Mr. 



is the photon energy in the Q rest frame and 



H 



^)(a? + a^) + (^^-^^ ' 



Pn -ps 

Ml 



)(ai6i + 0262 



2M= 



Mn' ' Mn 
Using the values obtained in the last section we get: 

T{n- S"7) = 5.6 X lO^^^eV 
which results in a branching ratio 

BR{n- S"7) = 7.0 X 10"^ 



^2) + -TT^('^l^l ~ '^262) 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



As mentioned before, the penguin contributions to this decay were calculated using 
bag model matrix elements, for weak and electromagnetic transitions at baryon level, and 
standard penguin coefficients, yielding p[: 



10 



gj^Pengmn^^- ^ ^-^^ _ 2 X lQ-^{Cp/Cif ~ 5 X 10-^ (4.6) 

where the ratio of the Cp/Ci coefficients of the QCD-corrected nonleptonic Hamihonian is 
approximately 1/20 [p] , ^ ,p^ . 

Our result suggests that the different contributions to the branching ratio BR{Q'^ — *■ 
S~7) arising from the single-quark s ^ dj transition, Eq. (|4.5|) , the penguin diagram, 
Eq.( [4.6| ), and the long distance process, Eq.( |1.4| ), are all comparable. In order to separate 
the contributions from these different processes it may be necessary to study the asymme- 



try parameters of the decay. It is also worth mentioning that the result given in Ref. HlO 



Eq.( |1.5|) , would be of the same order as our result, Eq. ( [4.5|) , had they used the same value 



of Crlrris) given here. In Ref. they also estimated another class of long distance con- 



tributions to the S~7 decay, using the vector meson dominance approximation, and 

they found that it could possibly saturate the experimental upper bound. However, their 
long distance result is proportional to a quantum chromodynamics coefficient, whose direct 
estimate is not reliable since all the calculation is far from the perturbative regime. For this 
reason we chose to quote here only the more traditional long distance result given in Ref. 

i- 

Our result is roughly a factor of 60 below the experimental upper limit, Eq. ( |1 . 1| ) , and it 
would be extremely interesting if new experiments could bring the upper limit down, since 
even adding all the contributions to this decay, one still gets a number that is at least one 
order of magnitude smaller than the experimental upper limit. 
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TABLES 

TABLE L Values of F2 and |Ai|F2 for different internal quarks. 
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FIGURES 



FIG. 1. Typical diagrams contributing to hyperon radiative decays at quark level. 

FIG. 2. Diagrams considered for the calculation of the Wilson coefficients of the correlator 
function. The crossed circle indicates the weak transition, with photon emission, and the cross on 
the s-quark line indicates that the correction to the propagator is relevant. 

FIG. 3. A qq pair in the vacuum coupled to the electromagnetic current at large distances. 

FIG. 4. The amplitude oi as a function of Mf for M| = (full hne), M| = 2.4GeV2 (dashed 
line) and M| = S.2GeV'^ (dotted line). 

FIG. 5. The amplitude ai as a function of M| for Mf = M| (fuh line),M| = lAGeV^ (dashed 
line) and M| = 2.0GeV2 (dotted line). 

FIG. 6. Same as Figure 4 for bi. 

FIG. 7. Same as Figure 5 for bi. 
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